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THC  COAGUI-ATIOH  OF  ClIAk'GED  CLOUD -DROPLIirrS 


L-  M.  Tx-'/in 

(Presented  by  0.  Yu.  Shmidt, 

Member  of  the  Acndeniy,  November  26,  1956.) 


According  to  present  theory  of  t,he  forination  of  precipitations 
[1,  2],  the  growth  of  cloud-droplets  in  the  range  of  dian^ters  d  =  4  ~  30  p 
takes  place  only  by  way  of  C(:)ndensation  of  water  vapor.  The  basis  of  this 
cGn-lusion  is  that  since  there  exists  a  critical  value  kcr  =  1.214  [2,  3] 
for  the  parameter  of  gravitational  coagulation,  then  coagulation  of  /two7 
neutral  droplets  is  jxissible  only  when  the  diameter  of  the  larger  drop  ~ 
exceeds  30  p.*  But  it  seems  to  us  that  growth  of  cloud-droplets  in  the 
range  d  =  4  ~  30  p  may  take  place  thrott5h  coagulation  of  electrically  charged 
droplets.  To  confirm  this  hypothesis,  let  us  cedculate  E,  the  coefficient  of 
captiare  (3j,  for  electrically  charged  droplets. 

Langmuir  [2]  and  Shishkin  [l]  unjustifiably  applied  the  theory  of 
deposition  of  an  aerosol  on  obstacles  to  the  theory  of  cloud-droplet  growth. 
The  coefficient  of  capture,  in  the  deposition  of  an  aerosol  on  u  sphere,  was 
calculated  [2]  under  the  assumption  that  the  sphere  is  considerably  larger 
than  tfie  aerosol  particle.  Consequently  these  calculations  may  be  applied 
to  the  descent  of  droplets  of  raindrop  size  in  a  cloud,  but  not  to  the  ca.se 
of  two  droplets  of  comparable  size.  Langmuir,  in  calculating  the  growth  of 
cloud -droplets,  really  considers  a.ll  the  cloud -droplets  to  be  fixed  in  sptce, 
except  one  droplet  (soraetlmeB  of  the  same  dimensions  as  the  others)  which 
moves  under  the  Influence  of  gravitation  and  of  rising  air  currents  [2). 
Ghishkin,  however,  takes  the  coagulation  of  droplets  as  depending  only  on 
their  relative  velocity  [1],  although  it  is  easy  to  see  that  the  air  veicclty 
fields  around  the  droplets,  which  determine  the  droplet  interaction,  depend 
on  the  velocities  of  the  droplets  relative  to  the  air.  Since  it  is  not  at 
present  possible  to  account  for  the  interaction  of  the  aerodynamic  fields  of 
two  droplets,  it  seems  reasonable  to  us  to  perform  the  calculations  for  cases 
when  one  of  the  drops  is  considerably  smaller  than  the  other,  and  In  doing  so 
to  take  the  serodynarolc  field  of  the  small  droplet  as  not  affecting  the  field 
of  the  large  droplet.  Then  in  a  coordinate  system  tied  to  the  large  droplet, 
the  small  droplet  will  be  acted  upon  by  the  force  of  gravity,  by  the  electric 
forces,  end  by  the  aerodynamic  Stokes  force. 


*  For  instance,  with  a  rising  air  velocity  of  10  cm/sec,  the  droplet 
diameter  grows  by  condensation  from  11  to  30  u  in  7.5  hours,  and  then 
grows  by  gravitational  coagulation  to  1300  u  in  15-20  minutes  (1). 
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l)  III  the  cane  of  both  droplets'  bearing  electric  charges  of 
different  signs,  the  eQoation  of  motion  of  the  small  droplet  will  take 
the  form 


jtrf*  d  V,  (r.) 
-5-P»-U7" 


•3r(irflu,  (ri)-  Vj  (I'l)]  -{- 


?,  ,  nrf* 

— J  r,  +  — Pkg, 
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where  r,  is  the  radius  vector  to  the  center  of  the  small  droplet ,  v,  =  c/fi/rf/ 
is  its  velocity  vector,  d  is  the  drojdet  diameter,  pjj  its  density,  ej,  eg 
the  charges  on  the  droplets*  u,  (r,)  is  the  field  of  velocities  of  the  air 
stream  around  the  large  droplet,  u  is  the  viscosity  coefficient  of  the  air, 
and  g  is  the  vector  acceleration  of  gravity. 


Converting  equation  (1)  to  the  diinenBionless  form  [3]  (as  characteristic 
velocity  we  take  the  constant  descent  velocity  of  the  large  droplet, 

“r,  —  and  as  cnaracteriEtio  "ize  D/2),  we  find  * 

-|-Vr=  U  +  feO,  ~  +g,,  (2) 

where 
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are  dimensionless  scaling  criteria  for  the  phenomena. 


It  is  obvious  that  for  |ai{  »  1  we  may  neglect  the  inertia  term  in 
equation  (2),  whereupon  this  equation  takes  the  fonn  (for  a  =»  -kcri  >  0): 


V==U-a'^+g,  . 


(3) 


Examination  of  equation  (3)  shows  that  for  small  values  of  a  it  has 
seven  singular  points  (four  saddles  on  the  p  axis,  two  centers  and  a  dipole 
at  tne  coordinate  origin),  and  for  large  values  of  5,  three  singular  points 
(two  saddles  on  the  o  axis  end  a  dipole).  Figure  1  illustrates  the  topo¬ 
logical  configurations  of  the  trajectories  of  equation  (3)  for  different 
valiies  of  a.  This  figure  shows  that  the  separatrices  A  pas  sing  through  a 
saddle  (p  •  0,  x  >  1)  represent  limiting  trajectories  determining  the  capture 
coefficient  E.  This  means  we  have  E  «  pi5  (Po  being  the  distance  of  the  said 
separatrlx  from  the  p  axis  at  x  -♦-«),  since,  as  we  shall  show  infra,  these 
oepara'crices  do  not  Intersect  the  surface  of  the  large  droplet. 


On  the  other  hand,  we  can  derive  expressions  for  the  trajectories 
described  by  equation  (3)  in  finite  form,  if  we  note  that  the  field 
corresponds  to  a  Stokes  stream  having,  in  dimensionless  (spherical)  coordinates, 

the  stream  function  ^  (r,fi)  Then  for  the  trajectories  we  get 


ir 


1  a 

- - 
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*  In  equation  (2)  we  have  the  dimensloalesB  quantities  v-v, 
r  -  2ri/;j,  ■r-.2/u^/.0  Where  D  is  the  diameter  of  the  large  droplet. 


tl  - .  u,  /  u„, 


and  consequently, 


jj,  j  { .  r*  sin*  k;\  * 

rfl-: -=  d  j — ^j==O5ln0d».  ^5j 

From  this  it  is  easy  to  find  the  equatir-n  of  the  trajectories: 

'h  ■=• 't'l.  »  —  0  {J  -f  cosO). 

where  is  the  value  of  for  the  trajectory  at  0  i»  «. 

Noting  that  the  separatvix  we  are  seeking  has  ,  at  0  •»  0,  a  finite 
value  of  r,  we  find  for  it  bhe  value  of  equal  .o 'l>i.  n  =  -  ■=2a.  Since 

according  to  equation  (5)  we  have  d^i/d9  >  0  along  U\e  trajectoiy,  then  the 
sep-'iratrix  in  question,  having  -  0  at  0  «  0,  nowh  re  intersects  the  surface 
of  the  large  droplet,  for  which  <  0.  Consequent]/, 


^  = 


*  “  gx  ' 


(7) 
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2)  Let  us  now  consider  that  case  when  the  large  droplet  Is  charged 
and  the  small  droplets  neutral.  Then  in  the  first  approximation  the  charge 
of  the  large  droplet  will  induce  dipoles  in  the  small  droplets  and  attract 
them.  The  equation  of  motion  of  the  small  droplets,  In  the  same  coordinate 
system,  takes  the  form  (ref.  [i]): 


<<K*  t/V,  (f,) 

"TP*  ~'dr~ 


—  3«iid  lu,  (r,)- V,  (r,)| 
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where  €  is  the  small  droplet  dielectric  constant. 


In  dimenslonleea  form  equation  (8)  becomes: 


*  +  (9) 


where  a. 


w  e  +  2 


>0. 


If  as  »  1,  then  equation  (9)  is  simplified; 

v«u  — — +  gi.  (10) 

and  the  trajectories  of  motion  may  be  determined  from  an  equation  similar 
to  (5): 


no,  sin  6  ltd 

- p - 


(11) 
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Equation  (10)  has  eight  singular  points.  Their  positions  and  the 
topological  configurations  of  the  trajectories  are  similar  to  the  cases 
illustrated  In  Figure  1.* 


-  oOo  - 


3)  Tables  1  end  2  present  the  results  of  our  calculations  of  E; 
likewise  values  of  CTj  and  03,  enabi  ixig  one  to  estiroati.;  the  applicability  of 
the  computation  methods  as  developed  (the  calculations  were  carried  out  for 
>  5  and  gx  ^  0.25).  The  order  of  magnitude  of  the  charge  was  cal¬ 
culated  according  to  Yn.  I.  Frenkel’  [S].,  and  for  the  abpolute  magnitvide  of 
the  charge  in  the  CGSE  system  we  accepted  the  expression 


e-0.5  lCr»rf(M,  (12) 

corresponding  to  a  value  of  0,3  V  for  the  electrokinetic  potential.  It  seems 
to  ub  that  the  magnitude  of  the  charge  accepted  by  us  according  to  expression 
(12)  has  been  selected  reasonably  as  regards  order  of  magnitude,  and  does  not 
contradict  the  experimental  findings  of  Gunn  [6),  who  measured  the  gross 
charge  of  cloud-droplets  of  diameter  d  >  10  p. 

TABLE  1 


OM 

d  U»1 

.  1 

■ 

•* 

5 

1 

15000 

65.2 

15 

6 

5.1 

2,8 

5 

2 

3900 

71.7 

20 

1 

61 

0.98 

7.5 

1 

3080 

18,8 

20 

3 

6,8 

1,00 

7,5 

*\ 

192 

25.9 

25 

1 

25 

0,50 

10 

1 

975 

7.9 

25 

2 

6,2 

0,50 

10 

5 

39 

10,1 

30 

1 

12 

0,29 

15 

1 

192 

2.3 

30 

1.5 

5,1 

0,29 

TABLE  2 


DM 

/,-o.s»| 

/r.-o.i 

|r.-o.o* 

e 

5 

1250 

2,32 

1.48 

0,89 

7,5 

110 

1.06 

0.69 

0.10 

10 

20 

0,57 

0.38 

0,22 

12,5 

5.1 

0,34 

0,23 

0,13 

*  The  difference  amounts  to  the  fact  that  the  dipole  reduces  to  two  nodes, 
a  stable  node,  at  the  coordinate  origin,  and  an  unstable  node  (p  0, 

0  <  X  <  1),  while  the  centers  are  converted  into  unstable  foci  or 
unstable  nodes.  As  previously,  the  separatrix  A  passing  through  a  saddle 
(p  =  0,  X  >  1)  does  not  intersect  the  sxirface  of  the  large  droplet 
(equation  (11))  and  constitutes  a  limit  trajectory.  Numerical  calculations 
have  shown  that  equation  (11)  is  convenient  for  constructing  this  separa¬ 
trix  (calculation  of  E). 


Tables  I  and  2  show  that  under  these  aBeumptions  the  coefficient  of 
capture  for  clcud-dropaets  ^  in  the  ran^ie  with  which  we  are  dealing,  is  high 
enough  to  guarantee  a  considerable  coaf^ulatior  of  droplets,  while  for  neutral 
droplets  of  the  same  "i^es  there  is  no  coagulation.  The  amount  of  coagulation 
will  be  rather  great  even  if  the  actual  size  of  the  charges  is  an  order  of 
magnitude  lower  than  that  which  we  have  taken. 


oOo  - 


4)  From  the  expressions  for  and  Oa  it  follows  that  at  large  values 
of  D  (when  u,„  also  is  large)  the  quantities  ctt.,  «  1  and,  consequently,  the 
effect  of  the  electric  charges  on  coag-olation  may  be  neglected.  This  acoo^ints 
for  the  experimental  results  of  Gunn  and  Hitschfeld  [7],  in  which  D  =  3.2  mm, 
d  >  5  p  and  *  0.2  CGSE.  Thus  in  these  experiments  u<o  »  825  cm/sec  and 
CTi  =  4- 10“^,  <  7*10"^.  Consequently  these  authors  did  not  find  any 

electric  charge  effect  on  the  coagulation  of  a  large  droplet  with  small 
droplets. 

In  the  same  way,  by  calculating  the  quantities  ai  and  Cz  it  may  be 
shown  that  in  cloud-di’oplet  catchers  the  precipitation  of  droplets  onto  the 
receiver  is  in  practice  independent  of  the  electric  charges  of  the  droplets 
and  of  the  droplet  receivei . 

Similar  results  are  indeed  found  for  charged  droplets. 

In  conclusion  I  express  icy  oo>rdial  thanks  to  V.  M.  Bovsheverov  for 
fruitful  discussion  cf  the  results  oi  this  work. 
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Academy  of  Sciences  of  the  USSR. 


REFEROCES 

1)  N.S.  SHISHKIN,  Usp.  Piz.  Rauk  ^^Progress  of  the  Riysical 

Sciencesy,  45,  3,  313  (1951). 

2)  I.  Langmuir,  J.  Meteor.,  5,  175  (1948). 

3)  L.,M.  LEV'IN,  ok.  Akad.  Nauk,  n,  No.  6  (1952). 

4)  R.  Cochet,  Ain.  de  G^ophys.,  8,  33  (1952), 

5)  Ya.  I.  FRENKEL',  Teoriya  yavlenii  atmosfemogo  elektrichestva 

^heory  of  atmospheric  electric  phenomenay,  1949. 

6)  R.  Gunn,  J.  Meteor.,  9,  397  (1952). 

7)  K.  Gunn,  W.  Hitschfeld,  ibid.,  6,  7  (1951). 


Translated  from  Dok.  Akctd.  Nauk,  ^  (1954);  G,  1045-1048. 


SIZE  DISTRIEirriON  FUNCTION  FOR  CLOUD-DROPLEl'S  AND  RAINDROPS 

L.  M.  Levin 

(Presented  by  G.  A.  Hajobur^ev, 

Kernber  of  the  Academy,  December  4,  1953.) 


Modern  theories  of  the  formatio:*.  rind  development  of  clouds  and 
forms  of  precipitation  az'e  essentially  t  .ncemed  with  empirical  sice  dis¬ 
tribution  functions  foi’  cloud -droplets  raindrops.  In  chronological 
order,  four  distribution  furictions  have  been  used  in  these  theories:  the 
Smolvkhovski -Schumann  theory  [1-3],  the  Marshall-Palmer  theory  [4],  the 
Best  theory  [5,  6]  and  the  Khrgian-Mazln  theory  [7].  Th(  first  of  these, 
a  theoretical  fornaila,  was  derived  by  omolukhovski  for  colloids,  and  the 
same  scheme  for  deriving  the  distribution  function  was,  without  any  suf¬ 
ficient  justification,  carried  over  by  Schumann  to  the  cloud-droplet  case. 
The  other  distribution  functions  ai'e  empirical. 

In  1941,  A.N.  Kolmogorov  proved  [3]  that  if  we  take  a  rather 
general  scheme  of  random  break-up  of  particles,  then  in  the  limit  we 
obtain  a  logarithmically  normal  size-distribution  of  particles.  It  is 
natural  to  suppose  that  a  logarithmically  normal  distribution  might  like¬ 
wise  be  found  for  the  random  coagiolntion  of  cloud-size  droplets,  in  the 
D.iinlting  case  of  a  great  number  of  particle -merger  events. 

In  1952-1953,  a  large  amount  of  steti'^ical  data  was  collected 
by  the  Elbruz  expedition  of  the  Geophysical  Institute,  Academy  of  Sciences 
of  T,be  USSR,  data  which  made  it  possible  to  examine  the  question  of  cloud- 
droplet  distribution  function  and  to  stats  the  premises  for  the  raindrop 
distribution  function^* 


Figure  la  shows  "rectified  diagrams”  [10]  for  the  size -distribution 
of  cloud- droplets  collected  from  Indivldtial  clouds.  Plotted  as  abscissae 
are  vhe  logarithnm  of  droplet  diameter  d  in  lal crons;  as  ordinates,  values 
of  [F(d)],  where  F(d)  is  the  selection  I'strlbution  function,  and 

(''^)  =  ^  {herd  is  the  inveesw'  rUiction  of  Y).  According  to 

•-.or. 

the  logarithmic  all;/  normal  distribution  Tnr  -.hich  the  density  of  distribution 
will  be  given  by  the  formula: 


rt{x) 


’V2ftx 


90* 


(1) 


The  experiir.ental  procedure  followed  by  the  expedition  in  its  reseeurch  on 
the  microstruc ture  of  clouds  is  described  in  reference  [9].  Rain  spectra 
were  studied  by  a  dyed  filter  paper  method. 


the  exr>_'rii:jc))iai  points  mst,  in  these  coordinates,  lie  on  a  straight  line. 

The  pc.int  cf  intersection  of  this  straight  line  and  the  axis  of  the  abscissae 
xs  In  and  tne  tangent  of  the  straight  line's  angle  of  slope  is  l/a  [10]. 

Figure  le  shows  tliat  the  ext^>eriinental  findings  for  clouds  are  satis¬ 
factorily  described  by  a  logarithmically  normal  distribution  function.  At 
the  present  time,  the  expedition  has  acc-nnulated  over  thirty  samplings  of 
cloud-droplets,  with  the  numbers  of  droplets  in  the  sample  ranging  from  2000 
to  20,000.  They  all  support  the  above  conclusion. 

is  for  the  raindrop  size  distribution,  we  are  not  at  the  present 
tipje  in  possession  of  any  similar  material  collected  from  individual  rain¬ 
falls.  Fiom  the  material  for  1952  it  was  possible  to  determine  only  the 
average  spectra  for  a  number  of  rainfalls  of  one  particular  intensity. 

These  spectra  are  shown,  ir,  the  same  coordinates,  in  Figure  2.  The  graphs 
of  Figure  2  indicate  that  rain-droplets  likewise  have  a  logarithmically 
normal  distribution  lav/.  Since  this  last  conclusion  is  dravm  only  for  the 
average  spectra,  It  is  of  great  importance  to  collect  statistical  material 
on  the  spectra  of  individual  rainfadls. 

Moreover,  to  answer  precisely  the  question  of  the  distribution 
fiuiction  It  is  essential  to  expand  the  range  of  captured  droplets,  both  for 
clouds  and  for  rain.  In  the  case  of  clouds,  it  would  appear  essential  to 
establish  a  procedure  for  capturing  fine  droplets  (d  <  4  u)  and  a  procedure 
for  collecting  a  greater  amount  of  the  large  droplets  (d  >  50-50  p).  In  the 
case  of  rain,  it  would  be  desirable  to  effect  the  capture  of  drops  of  dia¬ 
meter  d  >  400  w. 

It  seems  to  us  that  the  logaiithmically  normal  law  of  size-dist'^ibution 
of  cloud-droplets  In  the  diameter-range  d  ■  4  ~  50  u  proves  that  cr  rg-aiatlon 
of  cloud-droplets  occurs  over  the  whole  of  this  range.  This  fact  indirectly 
confirms  our  working  hypothesis  of  coag’.ilatlon  of  electrically  charged  cloud- 
droplets  In  the  dj ometer-range  d  »  4  -  50  p,  as  set  forth  In  reference  [11]. 

In  reference  [7]  it  was  shown  that  the  size-distribution  of  cloud- 
droplets  is  well  approximated  by  the  formula 

a(x)  -  ax^  (2) 


This  formula  was  derived  from  the  fact  that  over  a  large  range  of 


dlafieters  the  empirical  quantity  Ig 


Hill 

42 


was  found  to  be  linearly  related  to 


d.  But  this  relationship  holds  with  fair  accuracy  for  a  great  range  of  x, 
in  the  case  of  a  logarithmically  normal  distribution. 


In  Figure  3a  we  show  graphs  of  Ig  where  n(x)  is  given  by 

formula  (1).*  These  graphs  show  that  for  the  range  a  ■  0.25  -  0.5  (that  is, 
n  range  practically  covering  all  the  values  of  cr  encountered  in  tixe  clouds 
investigated),  relationship  (2)  Is  fairly  well  satisfied  for  x  >  xi  «  0.5 
1.0  X,,.  In  the  region  x  <  xi  there  is  a  characteristic  bend  in  the  curve  of 


The  quantity  Xq  was 
We  note  that  Xq  is 


taken  for  the  scale -unit  on  the 
the  median  of  the  distribution. 


of  Figure 
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'li — : 

'\Ad!Sf:j  S, 

4..; 

0, 

■  J 
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1  -  Sept.  3,  1953,  at  18:00  hre : 

N  =  27,800,  a  -  0.35,  do  -  8.6  n- 

2  -  Sept.  3,  1953,  at  18:30  hre: 

N  -  16,000,  cr  •  0.35,  do  -  9.3  n. 
•5  ••  Sept.  3,  1953,  at  19:00  hrs; 

N  «  7,500,  (T  -  0.31,  do  »  10.7  \i. 

4  -  Sept.  19,  1953,  at  10:45  hrs: 

N  -  4,352,  (T  -  0.30,  do  •  8.0  jl. 

5  -  Sept.  19,  1953,  at  19:00  hrs; 

N  -  3,900,  a  •  0.27,  do  «  0.8  |i. 


Pig.  2. 

1  -  I  •  5. 6-7. 5  saa/hr,  K  ■  28,700, 

<T  ■  0.52,  do  •  0.57  BSB. 

2  -  I  -  3. 6-5. 5  tm/hTf  N  «  57,900, 

0"  «  0.46,  do  •  0.64  BOB. 

3  -  I  ^  2. 6-5. 5  tm/hr,  H  ■  27,700, 

cr«i  0.48,  do  •  0.56  isn. 

4  -  I  -  1.6-2. 5  ncs/hr,  R  •  24,000, 

«r  -  0.46,  do  -  0.56  ssa. 

5  -  I  «  0.5-1. 5  im/hr,  N  •  22,300, 

or  ■  0.41,  do  “  0.59  Bsx. 

6  -  I  «  0.5-7. 5  mm/Vir,  N  ■  140,600 

<7  «  0.48,  do  «  0.58  nn. 


n(x) 


an 


j.g  • 


n(d) 

d2 


inflexion  which  is  seen  in  nearly  all  the  experimental  ciirves 
(see  Figure  lb)  corresponding  to  the  distributiona  shown  in 


Figure  la. 


Similaiily  it  may  be  shown  that  the  size-distribution  of  raindrops 
as  found  empirically  by  Marshall  and  Palmer  [4]: 

n(x,)  <»  ae  (3) 

is,  over  a  large  part  of  the  curve,  an  approximation  to  the  logarithmically 
normal  distribution.  This  fact  is  illustrated  in  Figure  3b,  which  is  a 
graph  of  Ig  n(x)  for  n(x)  as  given  by  formula  (1). 

Here,  incidentally,  we  must  note  that  the  experimental  data  obtained 
by  us  on  cloud-droplets  do  not  satisfy  the  Saolu!<hov8ki-Schumann  distribution, 
for  which 

n(x)  *  ax^  (4) 

This  follows  from  the  fact  that  in  the  coordinates  of  Figixre  lb  this 
distribution  (4)  should  appear  as  a  cubic  parabola.  As  may  be  seen  from 
Figure  lb  this  is  not  the  case. 

In  conclusion  we  express  our  cordial  thanks  to  A.M.  Obukhov, 
Member-Correspjondent  of  the  Acadeny  of  Sciences  of  'the  USSR,  for  his  valuable 
advice.  > 
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